Spontaneous symmetry breaking in a bridge model fed 
by junctions 



Vladislav Popkov 12 , Martin R. Evans 3 and David Mukamel 4 

1 Dipartimento di Fisica "E.R. Caianiello", and Consorzio Nazionale Interuniversitario per 
le Scienze Fisiche della Materia (CNISM), Universita di Salerno, Baronissi, Italy 

2 Interdisziplinares Zentrum fur Komplexe Systeme, Romerstrasse 164, D-53117 Bonn, 
Germany 

3 SUPA, School of Physics, University of Edinburgh, Mayfield Road, Edinburgh EH9 3JZ, 
UK 

4 Department of Physics of Complex Systems, The Weizmann Institute of Science, 
Rehovot 76100, Israel 

E-mail: popkovOsa. inf n . it ,m. evansOed. ac . uk,david.mukamel@weizmann. ac.il 

July 28, 2008 

Abstract. We introduce a class of ID models mimicking a single-lane bridge with two 
junctions and two particle species driven in opposite directions. The model exhibits 
spontaneous symmetry breaking (SSB) for a range of injection/extraction rates. In this 
phase the steady state currents of the two species are not equal. Moreover there is a co- 
existence region in which the symmetry broken phase co-exists with a symmetric phase. 
Along a path in which the extraction rate is varied, keeping the injection rate fixed and 
large, hysteresis takes place. The mean field phase diagram is calculated and supporting 
Monte-Carlo simulations are presented. One of the transition lines exhibits a kink, a feature 
which cannot exist in transition lines of equilibrium phase transitions. 
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1. Introduction and Model Definition 

Nonequilibrium stationary states (NESS), in which probability currents are supported 
and invariant measures are not generally of Gibbs-Boltzmann form, offer many surprising 
phenomena not seen in equilibrium systems. Examples include boundary-induced phase 
transitions and spontaneous symmetry breaking (SSB) in one-dimensional systems. The 
minimal models for NESS are driven particle models such as the Totally Asymmetric 
Exclusion Process (TASEP) and its relatives. Some exact solutions for these models have 
been found, for example the TASEP, leading to an understanding of boundary-induced phase 
transitions pp. However, important questions remain as to the nature of SSB in driven one- 
dimensional systems. 

The SSB phenomenon was first observed in what is often referred to as the bridge 
model [2], and subsequently in some other models [3j H]. The bridge model comprises a 
one dimensional lattice on which positive particles (pluses) move to the right and negative 
particles (minuses) move to the left. At the left boundary of the lattice, pluses may enter 
the lattice and minuses may leave; at the right boundary of the lattice minuses may enter 
the lattice and pluses may leave. When the exit rate at which particles may leave the lattice 
is lowered, the stationary state changes from a symmetric one, in which the currents and 
bulk densities of pluses and minuses are equal, to a symmetry-broken state where there is 
a majority species of particle with larger current and higher bulk density than the minority 
species. However, the symmetry broken state has not been solved exactly except in the limit 
in which the exit rate tends to zero where it has been rigorously proven that SSB occurs[5]. 

The transition to this symmetry broken state remains a subject of debate. A mean 
field theory originally predicted that the transition should occur via an intermediate, weakly 
symmetry-broken phase j6]. Monte Carlo simulations have shown that the mean field theory 
does not correctly predict the position of the transition [7] , however, at least on finite systems, 
an intermediate phase is seen and the transition from symmetric to the strongly symmetry- 
broken phase occurs through a sequence of transitions [8]. This sequence has also been 
observed in a related model [I] . However, it has been suggested that the region of parameter 
space occupied by the intermediate phase and over which the sequence of transitions occurs 
disappears in the infinite system size limit [9]. 

The failure of the mean-field approximation to exactly predict the phase diagram can 
be traced back to the boundary conditions of the bridge model which do not correspond to 
particle reservoirs at fixed density. Instead there are effective impurities at the boundaries 
[10]. Also the hydrodynamic limit of the bridge model is not well-defined over the whole 
phase space. In particular, it breaks down at the SSB transition in the limit of small input 
and exit rates [10]. 

In this paper we introduce a new class of bridge models demonstrating SSB phenomenon. 
In these models the input and output of the pluses and minuses are governed by TASEPs. 
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Thus one can think of the ends of the bridge as junctions where TASEPs for the pluses and 
minuses merge. The input and output rates at the bridge are not external parameters as 
in TASEP, rather they are determined self consistently by the dynamics of the bridge and 
its feeding segments. The models exhibit complex phase diagrams including a phase co- 
existence region showing hysteresis and quite special triple phase co-existence point, which 
appear to be correctly predicted by a mean field theory. 

Our model is defined as follows. We consider a chain of the length L occupied by two 
species of particles, "plus" particles moving to the right and "minus" particles moving to 
the left, with hardcore exclusion and random sequential update. The chain is region II of 
Fig. [T]and we shall refer to it as the bridge. At each end of the bridge, there are junctions 
where the chain splits into two parallel segments, one containing only plus particles and holes 
and the other only minus particles and holes (sections I and III in FigJT]). At the external 
boundaries of the parallel segments, usual TASEP rules are applied. Namely, plus particles 
are injected into their segment of section I at the left with rate a, if the first site is empty, 
and are removed with rate (3 from the right end (upper segment of section III in Fig. [1]), if 
the last site is occupied. Likewise minus particles are injected into their segment of section 
III with rate a and removed with rate (3 from the left end of their segment of section I. In 
sections I, III pluses hop to the right inside their segments with rate 1. Plus particles enter 
the bridge at the left junction if the first site of section II is empty and leave it at the other 
junction if the first site of the plus segment of section III is empty, both with rate K. Inside 
the bridge, particles exchange with empty sites and with the minuses with the same rate 1. 
Similarly, inside the lattice minus particles hop to the left with rate 1 in sections I and III, 
enter and leave the bridge if the entrance or exit site is empty with rate K, and inside the 
bridge exchange with both pluses and empty sites with rate 1. The model is symmetric with 
respect to simultaneous charge inversion and left-right reflection. 

In this work we consider only K — 1, leaving a and (3 as the model parameters. We 
calculate the phase diagram of the model using a mean field approximation and direct 
simulation of the dynamics. 

2. Phase Diagram and Phase Transition Lines 

Before discussing the phase diagram of our model, it is useful to recall the phase diagram 
of the TASEP (for plus particles moving to the right). The phases are distinguished by the 
large system size limits of the expressions for p the bulk density of particles far from the 
boundaries and j = p(l — p) the current of particles. When a < 1/2 and (3 > a the low 
density (LD) phase occurs where the density in the bulk is equal to a and is determined by 
the left boundary; when (3 < 1/2 and a > (3 the high density (HD) phase occurs where the 
density in the bulk is equal to 1 — (3 and is determined by the right boundary; when a > 1/2 
and (3 > 1/2 the maximal current (MC) phase occurs where the bulk density is 1/2. 
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Figure 1. The bridge model with two junctions. Positively (negatively) charged particles 
hop to the right (left). The model is invariant with respect to left-right reflection and 
charge inversion. Section II is the bridge. It contains positive and negative particles and 
holes. Sections I and III comprise parallel segments each containing pluses and holes or 
minuses and holes. 
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Figure 2. Mean field phase diagram of the bridge model fed by junctions for K = 1 (see 
text for the definition of the phases). The dotted lines show the location of the LDS1-LDS2 
and LDS2 — > SSB phase transitions from the Monte-Carlo simulations. 



We now present the phase diagram of the bridge model fed by junctions obtained from 
a mean field analysis detailed in the following section. The resulting a- (3 phase diagram is 
given in Fig. [2j it contains three different phases. 

(a) Low density symmetric (LDS1) phase (a < 1/3,/? > a). Each species establishes 
a homogeneous state with low particle density p = a in each segment, see FigJHJ^a). 

(b) Spontaneous symmetry-broken (SSB) phase (/? < 1/3, a > (3). The two species 
have different densities and fluxes and the phase comprises two symmetry related states. 



1 



The majority species (the pluses in Figj3]^b)) establishes a high density state with bulk 
density 1 — (3 in all segments while the minority species (the minuses in FigJ3](b)) has 
bulk density (3/2 in the bridge and section I, and bulk density 1 — (3/2 in section III. 
(c) Low density symmetric (LDS2) phase (a > 1/3,(3 > 1/3). The pluses have bulk 
density 2/3 in section I and bulk density 1/3 in sections II, III whereas the minuses have 
bulk density 1/3 in sections I, II and 2/3 in section III. Thus the profile of the minuses 
mirrors that of the pluses, see Figj3](c), and the phase is symmetric. Note that while 
both densities on the bridge (section II) are low, on the other sections I, III the density 
of one of the species is low while that of the other species is high. This phase is thus 
different from the LDS1 phase. 

The phase diagram (Fig. [2]) is reminiscent of the TASEP phase diagram discussed above, 
but with several important differences. First the LDS2 phase replaces the maximal current 
phase and the transition lines to the LDS2 phase are at a = 1/3, (3 = 1/3 rather than 
a = 1/2, (3 = 1/2. Secondly, the SSB phase replaces the high density phase of the TASEP. 
Clearly a high density symmetric phase is not possible in the present model since in that 
phase the densities of both species on the bridge section have to be greater than 1/2. Thirdly, 
there is a co-existence region (1/3 < (3 < 2/5, a > (3) where the system may be in either of 
the LDS2 and SSB phases; this produces an interesting kink in the phase boundary of the 
SSB phase. 

We now discuss the phase transitions and transition lines separating the phases described 
above. For the TASEP, as well as other nonequilibrium steady states with a non-zero con- 
served current, the 'order' of a phase transition is determined by the non-analyticity of the 
current at the transition, i.e. the order of the derivative of the current at which a singularity 
appears pQ. In the present model we find some novel features for the nonequilibrium phase 
transitions. Notably, across all the transition lines some of the bulk densities change discon- 
tinuously, as can be seen by comparing stationary densities in different phases. Indeed for 
the LDS1/SSB transition there is even a discontinuity in the current of the minority species. 
We now describe the behaviours at the transitions. 

LDS1/LDS2 transition line: a = 1/3, (3 > a. The bulk densities in the bridge are 
continuous across the transition and the currents are continuous. In section I the plus den- 
sity jumps disco ntinuously across the transition from 1/3 to 2/3, whereas the minus density 
is continuous. Similarly in section III, the minus density jumps discontinuously from 1/3 
to 2/3. On the transition line one finds a shock in section I separating regions of plus den- 
sities 1/3 and 2/3 and a shock in section III separating regions of minus densities 2/3 and 1/3. 

LDSl/SSB transition line: a = (3 < 1/3. The bulk densities in all sections are discontin- 
uous across the transition: the bulk density of the majority species jumps from a to 1 — a in 
all sections; the bulk density of the minority species (taken as minus) jumps from a to a/2 
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Figure 3. Average density profiles for pluses and minuses, from Monte-Carlo simulations, 
in the LDS1 phase (Panel (a)), in the SSB phase with pluses majority (Panel (b)) and in 
the LDS2 phase (Panel (c)). Pluses (minuses) correspond to closed (open) circles. The 
system of 300 sites was equilibrated and then, averaging over 10 6 Monte-Carlo steps was 
done. Parameters: (a) a = 0.2, j3 = 0.4; (b) a = 0.2, = 0.25; (c) a = 0.8, (3 = 0.9. 



in sections I and II and from a to 1 — a/2 in section III. The majority current is continuous 
whereas the minority current jumps from a(l — a) to a/2(l — a/2). On the transition line 
two types of shock configurations are in fact observed. First there is symmetric configuration 
where the plus density is a in sections I and II, and in section III there is a shock separating 
regions of plus densities a and 1 — a. Similarly the minus density is a in sections II and III, 
and in section I there is a shock separating regions of minus densities 1—a and a. Secondly, 
there is an asymmetric shock configuration where, taking the majority species to be plus, 
the bulk plus density is 1 — a in sections II and III and in section I there is a shock between 
regions of plus density a and 1 — a. The minority species density, taken to be minus, is a/2 
in sections I and II and 1 — a/2 in section III. Here the currents of plus and minus are unequal. 

LDS2/SSB Co-Existence Region: a > f3, 1/3 < /3 < 2/5 . In this region, denoted 
"SSB+LDS2" in the phase diagram Figj2], the system can be in either of the two symmetry- 
related SSB states or the LDS2 phase. In the stochastic model of infinite size, the three 
phases are stable. In finite systems flips between the phases take place, with typical time 
between the flips growing exponentially with the system size L (this issue will be addressed 
in SecJH). Note that spatial co-existence of the phases (like e.g. in [3]) is not possible since 
they carry different currents. 
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Figure 4. Hysteresis path for the majority density in the bridge section in the coexistence 
region: keeping the injection rate constant a > 2/5, the extraction rate j3 was changed from 
to 1 and back. The majority density in the bridge segment is shown as function of (3. On 
increasing (3, SSB phase is stable for (3 < 2/5, while on the way back (decreasing 0) the 
LDS2 phase (p = m = 1/3) is stable for 1/3 < j3 < 1. The (barely visible) dotted line is the 
result of Monte-Carlo simulations 



In the infinite system, the co-existence region SSB+LDS2 entails hysteresis, which 
is observed when keeping the injection rate a > 2/5 constant and increasing/decreasing 
adiabatically the extraction rate (3 from (SSB phase) to 1 (LDS2) and back. Indeed, along 
the path — ► 1 the SSB/LDS2 transition happens at (3 = 2/5, while on the backward path 
it will happen at j3 = 1/3. The hysteresis is illustrated in FigJU 

A kink point with coordinates a = /? = 2/5on the phase diagram is an unusual one. 
In phase diagrams of systems at equilibrium, such a kink would imply that there must be 
another first order transition line emerging from that point, making it a triple point. This 
is a result of the Clausius-Clapeyron relation. Our model is out of equilibrium, and the 
existence of a kink, although unusual, does not violate any rule. 

3. Mean field solution 

We now derive the phase diagram through a mean field approximation commonly used for 
one-dimensional stochastic systems, wherein two-point correlation functions are replaced by 
products of one-point correlation functions. For the TASEP this mean-field approximation 
is known to predict the correct phase diagram and stationary bulk densities [H] . 

Let us denote the stationary density of plus particles at site k as Pk- Then, the mean 
field approximation gives the stationary flux of pluses as 

j+ =p k (l -Pk+i) =PsO- ~Ps) (1) 
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where S'=I,II,III and p$ is the limiting bulk density far away from the boundaries in segment 
S. Similarly, the density of minus particles at site k is m k and we write 

j_ = (1 - m k )m k+1 = m s (l - m s ) . (2) 

In the stationary state the bulk densities ps, rris can be different in different segments S but 
the currents have the same values everywhere. At the boundaries the currents read 

3+ = <x(l - pi) = Ppsl (3) 
j_ =a(l - m 3L ) = (3mi (4) 

and at the junctions, 

3+ = Pl(1 ~ Pl+i - m L+ i) = p 2L (l - P2L+1) (5) 
3- = m 2 L+i(l - P2L - m 2L ) = m L+ i(l - m L ) . (6) 

Following P,d2], we define effective entrance rates oig , a.g and exit rates ft, ft for each 
segment. For plus particles 

+ + 3+ + 3+ /^\ 
aj = a; oT n = - ; aj n = = p 2L ; (7) 



Pl+i J- - P2L+1 

- m L+1 ; ft } = — 

PL P2L 



(3t = — = l-p L+ i-m L+1 ; ( 3+=^ = i- P2L+i; (3+ n = f3. (8) 



The rates a^,/?^ for minus particles are obtained by the substitutions / III and 
Pk — > msL-k+i- For each species, each segment S can be viewed as a TASEP model with, e.g. 
for the pluses, the effective injection/extraction rates a^, (3$. In the large segment length 
limit, the bulk density of pluses ps in segment S can be read from the phase diagram of a 
TASEP H3III3, namely 

LD phase < (3$ , < 1/2 ps = , (9) 

HD phase 1/2 p s = 1 - ft , (10) 

MC phase a^,ft>l/2 p s = 1/2 . (11) 

In the HD phase the density profile is flat at the right boundary and in the LD phase the 
density profile is flat at the left boundary. For minuses, the corresponding bulk densities 
are given by (T9 H1TT) with substitutions + — > — and p — > m. Thus, in the large segment 
length limit, each segment of the bridge model must exhibit one of the solutions (EHUJ). The 
possible solutions of the mean-field equations [TJ-[6] (MFE) are listed below. 

(a) LDS1 phase. Each segment is in the low density phase ([9]) with density p$ = ms = a 
and currents j + = j- = a(l — a). The LD density profile implies that Pl+i = P2L+1 = 
m 2L = m L = a and the symmetry of the phase implies that m i+1 = p 2 L, rn 2 L+i = Pl- As 
the solution is symmetric we need only consider the effective rates a g, ft: we read off from 
((?]) afj = afjj = a and from (jSJ) ft = 1 — Pl+i — rn L+1 = 1 — 2a, ftj = 1 — a. Conditions 
(Q for each segment reduce to the key conditions 

a<p, a< 1/3 (12) 
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The first inequality follows from af H < (3^ n and the second follows from af < (3f . 

(b) SSB phase. Let us consider the solution of the MFE where the pluses establish 
a high density phase pg = 1 — (3 > 1/2 in all segments. The minuses (the minority 
species), are in a low density phase m$ = 7 in segments S = I, II and a high density 
mm = 1 — 7 in segment III giving a current j_ = 7 (1 — 7) < j + . The value of 7 is to be 
determined. The structure of the profiles imply that pi = P2L = 1 — P, ™>l = m 2L = 7, 
m 2 L+i = 1-7- Then (06]) imply that p 2L+1 = 1-/3, m L+1 = 7, p L+1 = 1 - /3 - 7. /From 
one has j_ = 7 (1 — 7) = (1 — 7) (1 — (1 — (3) — 7) giving 7 = /3/2. We read off from ([7j) 
a/j = 2(1 — /3)/3, af n = 1 — /3 and from (jHJ) = = /3. Similarly we find «7 = ajj = /3/2 
and f3jj = 1— P/2, Pj n = (3/2. Conditions (jUj) for each segment reduce to the key conditions 



the first coming from Pf < a\ and the second from (3^ < afj. By a symmetry 
transformation (spatial inversion and interchange of minuses and pluses) one obtains the 
other symmetry-broken solution with minuses in the majority. 

(c) LDS2 phase. In this solution each species is in a low density phase in the bridge 
and the segment fed by the bridge and is in a high density phase in the segment feeding the 
bridge. Thus pi = 1 — 5, mm — 1 — 5, p$ — m$ = 5 (remaining segments) where the value 
of 5 < 1/2 is to be determined. As this solution is symmetric we need only consider the 
effective rates atg, (3$. The structure of the profiles implies thatpz, = 1 — 5, pi+i = P2L+1 = 5, 
m L = m 2L = 8, m 2L +i = 1-6. Then © gives j + = (1 - 5) (1 - 5 - p 2L ) = p 2L (l - 5), 
implying that p 2 L = 5 and 5 = 1/3. We read off from (JTj) a/ 7 = af n = 5 and from (jS|) 
(3f = 5, Pfj = 1 — 5. Conditions Qj for each segment reduce to the key conditions 



the first coming from (3j < aj and the second from aj n < (3j n . 

The question of the accuracy of the MF solution remains. Discrepancies between the 
predicted stationary bulk densities and Monte-Carlo simulations lie within numerical error 
bars, except for the LDS2 phase. In this phase the bulk density from Monte Carlo simulations 
Pn ~ 0.339 appears to differ slightly from the one predicted by mean field pn = 1/3. This 
discrepancy derives from the fact that while the LDS1 and SSB phases are controlled by 
uncorrelated inputs from the outer boundaries, the LDS2 phase is controlled by the injection 
rate from the junctions which could imply correlations in the input. It appears that such 
correlations result in only a small shift of the LDS2 «-»• LDS1 and LDS2 — > SSB phase 
boundaries to a ~ 0.339 < (3 and f3 ~ 0.339 < a, while other phase boundaries remain as 
in FigJ21 In view of this, we believe that the mean-field predictions are qualitatively correct 
and are quantitatively very accurate. 



a>/3, (3 < 2/5 , 



(13) 



a > 1/3 , p > 1/3 



(14) 
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Figure 5. (a) Average density of pluses in middle segment versus time, inside the SSB 
phase ( Panel(a)), and inside the SSB+LDS2 phase (Panel (b)). Panel (a): flips between two 
quasistable states p = 1— (3 andp = (3/2 are seen. Parameters are a — 0.45, j3 = 0.26, L = 20. 
Panel (b): Flips between three quasistable states p = 1 — /?, p = [3/2 and p = 1/3 are seen. 
Parameters are a = 0.45, (3 = 0.37, L = 100. 



4. Mechanisms of flips between symmetry-broken states in the SSB phase 

In a finite system, flips between the two symmetry-broken states (majority and minority 
species interchange) are observed (Figfj^a)). To demonstrate spontaneous symmetry 
breaking in the SSB phase and the SSB+LDS2 region FigJ2l we need to show that the flipping 
times between the stable states in finite systems diverge exponentially with the system size 
L. For the SSB phase the mechanism is best illustrated in the SSB region with the additional 
condition a < 1/2. In this region the plus profile in section I of Figj3]^b) is produced by 
a domain wall (a shock) between a region adjacent to the boundary with density a and a 
region with density 1—/3. The domain wall is biased to the left and therefore sticks to the left 
boundary. (Note that the motion of the plus domain wall does not affect the state of minus 
particles.) In the SSB configuration the entrance of minuses to the bridge is partially blocked 
by the high density 1 — (3 of pluses. To flip the SSB configuration Fig 13th), the domain wall 
of pluses pi — 1 — (3 sticking to the left boundary must retreat (against its bias) back to 
segment III in order to give an opportunity for the minuses to take over the key middle 
segment, block the entrance of pluses on it, and make the flip. The time to wait for such an 
improbable event grows exponentially with the system size L. Roughly, one can estimate the 
flipping time as follows. The shock of pluses (a, 1 — (3) is driven to the right with the rate 
r = (3(1 — (3)/(l — (3 — a) and to the left with the rate / = a(l — o)/(l — (3 — a). In the SSB 
region (fl~3|) / > r. Let us represent the shock position with a phantom particle which hops 
to the nearest left/right site with the rates I and r. Introducing the stationary probability 
rjk for the phantom particle to be at site k, where k — 0, 1, . . . 2L we find r\ k l = r] k _ir so that 
Vk = Vo ( r /l) k - We can estimate the flipping time by l/r] 2 L [IS], giving 

oor (r\- 2L / t \ i , a(l — a) 

tflip ~ (j j = exp (kL) where « = 2 log v . 
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In the rest of the SSB region (for a > 1/2) the mechanism is similar, but the density 
associated with the left hand boundary is 1/2 rather than a. Consequently, in the flipping 
time estimate a(l — a) has to be substituted with 1/4 for this case. 

For the region SSB+LDS2 with three stable phases, flips between the steady states 
are also caused by the shock motion, but their mechanism is different. Firstly, there 
is no direct transition between the two SSB states, only through the intermediate phase 
SSB1<^=^LDS2^=^SSB2, see Figfj^b). The existence of the intermediate LDS2 state makes 
the transition much easier: indeed, to flip from the LDS2 to SSB state, a shock need only 
cross one segment as opposed to two segments in direct transitions between the two symmetry 
broken states. This explains the big difference in system sizes between Figj5]^a) and Fig|5]^b). 

5. Conclusions 

In this work we have introduced a bridge model fed by junctions where the input and 
output streams of the bridge are themselves TASEPs. This has allowed the phase diagram 
predicted by the mean-field approximation to be confirmed by numerical simulations as 
being qualitatively correct and rather accurate quantitatively. The phase diagram exhibits a 
number of novel and interesting features that we discussed in Section II. We note here some 
further points that will be pursued in future work 

A particularly interesting feature of the phase diagram is that we have two types of first- 
order transitions: one with a co-existence line (a = (3 < 1/3) and one with a co-existence 
region (marked LDS2+SSB in Fig. [2]). In both cases the co-existence is between different 
possible steady states for the systems. The dynamics of how the system flips between these 
steady states is rather intricate but can be understood in terms of domain wall dynamics 
(work in progress). 

In the phase diagram Fig.[2]a phase resembling the maximal current phase of the TASEP 
is absent. We also note that none of the solutions of the MFE supports maximal current 
j + = j- = 1/4 in the system. The reason is that for the case K — 1 considered here, the 
junctions act as effective bottlenecks. For larger K (more precisely, for K > 2), a new phase 
reminiscent of the MC current phase in the TASEP does appear. 
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